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CHAPITRE III 



III APPLICATIONS DES EQUATIONS GENERALES A L 1 ETUDE 



D'ECOULEMENTS DE FLUIDES INCOMPRESSIBLES 



Pour resoudre les problemes pos£s, nous avons a notre disposition 
5 relations t 

la 

- Equation de continuity : traduit la conservation de'masse 

- Equations d'EULER : Equations dynamiques 

- equation d'Etat du fluide 

- theoreme des quantites de mouvement ou theoreme d'EULER 

- theoreme de BERNOULLI. 

Ces relations ne pourront etre appliquees a des phenomenes physi- 
ques que si les hypotheses relatives a. ces relations sont respec- 
tees. II est par consequent n^cessaire d* avoir une certaine con- 
naissance des phenomenes, connaissance qu'on peut acquerir par 
1' experience. La quality d'une theorie pourra etre testee par 
1* experience. 

III-I Applications du theoreme de BERNOULLI : 

III-I I) Pression statique. dynamique et totale s 

Reprenons 1* equation de BERNOULLI : 





Remarque 



On confond souvent p et p+^gz, les forces de pesanteur 
etant souvent negligeables par rapport aux forces de 
pression. 
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Tube ri^zometrique : 

Considerons une conduite assimilable h un tube de courant. Appli- 
quons la formule de nERNOULLI . 2 

P + K «h + - Gte - P t 

Entre 2 sections, on a done i 0 

V I V 2 
P I *t gh I — = P2 + f« h 2 + C"T~" 

Si la conduite est cylindrique ou 

prismatique, la vitesse est paralle- 

le aux generatrices et les lignes 

de courant sont des droites // aux 

generatrices. 

Dans ces conditions, la repartition des pressions suivant la nor- 
male aux lignes de courant est hydrostatique. Pour une section 
droite p = p + Pgh = Cte ■ pression statiqus 

La pression statique peut^tre mesur6e a l'aide d'un tube pl6zome - 
trique . C'est un tube de forme et d* inclinaison queloonquts qui 
debouche dans la conduite et qui est ouvert aux deux extremites. 
L'extremite donnant sur la conduite est une prise de presBion sta - 
tique . _ , « Le fluide monte dans le tube jus- 

qu'en Bp la cote de Bj etant £ga- 

En effet p+£gh est constant le long 
du tube^^hydrostatique) . Cette 
quantity est aussi constants le 
long de la section droite S de la 
. . , conduite. Si la pression ns subit 

pas de discontinuity au niveau du tube, p+^gh ne change pas en Aj 
quand on passe du tube dans la canalisation. 
On a done p + fg h = p„ + f g h 

Aj Aj flj V Bj 




B I 
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B I 



La cote de jB 



+ h , p etant la pression relative 



mesure done bien 



+ h regnant au niveau de la sec- 




tion S. Elle ne depend pas de la position de la prise de pression 
sur le pourtour ou a l'interieur de la section droite de la ca- 
nalisation. 

Remarque : I) 1 ' introduction de la prise de pression ne doit pas 
perturber 1 'ecoulement . Sur le pourtour, on peut uti- 
liser une petite cavit£ C qui set remplit 
de fluide mort. A l'interieur d'une cana- 
lisation, on associe la prise de pression 
statique a une portion de paroi plan conti- 
nue parallele aux lignes de courant. 
2) p est la pression relative : pression 
prise a partir de la pression atmospheri- 
que p = p' - p a 

III-I 2) Ecoulement a travers un orifice perce" dans une paroi 
mince - Theoreme de TORRICEILI ; 

Gonsid^rons un reservoir de section S comportant un orifice cir- 
culaire de section s. L'orifice est a bord mince t il est taill^ 
en biseau vers l'ext^rieur de sorte qu'on peut consid^rer qu'on 
a affaire a une paroi mince. 

Le liquide s'6coule h travers 
l'orifice et forme dans l'air un 
jet vertical. Par visualisation 
du liquide, on peut observer les 
lignes de courant 8 tout le li- 
quide participe a 1' ecoulement 
et les lignes de courant conver- 
gent avant d'atteindre l'orifice 
et cette convergence se poursuit 
au-dela de l'orifice. 
La section du jet diminue rapidement jusqu'a atteindre une section 
constante s_ qu'on appelle section contractee . Les lignes de cou- 
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rant k partir de s c sont pratiquement paralleles. 

Le rapport C 0 = s 0 / 8 * st ^ e coefficient de contraction . 

Ce coefficient est voisin de 0,6 pour les orifices circulaires ou 

memos rectangulaires. Si l'orifice a le profil d'une tuyere le jet 

libre ne pr£sente plus de contraction et 

dans ce cas s = s et C =1. 

c c 

Quand on se d^place perpendiculairement 
a des trajectoires reptilip-nes ou peu cour- 
tages (rayon de courbure R grand), la com- 
posante de 1 ' acceleration suivant la 
normale ~n princ^pale a la trajectoire est 
nulle ( * 0) 

Projetons 1' equation d'EULER suivant cette normale T ; il vient t 

projjj(F*- grad p) = 0 
Dans le champ de pesanteur "T = - gracl V = - gra3 gh 
soit projj^gracUp+Pgh) •» 0 
grt3p*\ dlT- 0 
ou - ^P =0, p* « Cte dans une direction normale aux 
trajectoires ou lignes de courant. 

La repartition des pressions dans une section normale aux lignes 
de courant quand celles-ci sont droites ou peu courbees est hydro - 
statjque . 

Cette condition n'est pas r^alis^e dans la partie convergente du 
jet. Par contre pour la surface s c , on a (p + £gz) c = Cte 
Mais pour la surface s c , z c = Cte done 

p m cte = pression au bord du jet dans la section s c = p fi 
Conslderons l'ecoulement pendant un temps assez court pour que la 
variation du niveau du reservoir puisse 6tre n^glig^e et que le 
mouvement soit consider comme permanent. Conside>ons un lilet 
fluide issu d'un point M de la surface libre et passant par un 
point C de s c { ngus pourrons applique jj le theoreme de BERNOULLI : 



y— + g + ^ 

s p 



Pc 

+ « Z C + ~f 
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I 

~T 



(v 2 c - v*) = 



+ « (z M " 



Z c> 



M -z c = h et si on neglige la vitesse a la surface li- 

Je debit en vo- 



Bi on pose z 

bre (etant donne* le rapport des sections S et s ( 
lume dq du filet de courant est 4gal a dq = V„ 

v v q vi 

dS„/ds est du mime ordre de grandeur que — — - et par consequent 
V M c s o 

), on obtient pour V c : 



dS M . V c ds c . 



2 ( 



"I 



+ gh) 



(I) 



Examinons un cas particulier important : Le recipient est en com- 
munication avec 1' atmosphere ; p M = p et la relation (I) devient : 
i « a 

V c =y2gh l (2) 
c'eat le th^orfeme de TORRICELLI ♦ 

La vitesse du fluide non visqueux dans la section contractee est 
la mime que s'il tombait en chute libre d'une hauteur h. 
REMARQUE S t I) La forme du jet apres la section contractee n'est 

pas exactement cylindrique mais legerement conique (cf 
exercice s application du thSoreme de BERNOULLI) 

2) Orifice a paroi mince perce sur la paroi verticale . 

La ligne moyenne du jet n'est plus 
verticale mais forme une parabole 
d'axe vertical dont la position per- 
met de determiner la vitesse dans la 
section contractee. Mais le jet pre- 
sente^une section contractee. 
Les resultfts precedents peuvent etre 

appliques a condition que le diametre 
de 1* orifice soit assez ftible pour que 1 ' on puisse admettre que 
z c = Cte. 

3) L' application du theoreme de BERNOULLI n'est valable que si le 
fluide n'est pas visqueux^ Toutefois s'il n'en est pas ainsi, le 
frottement du fluide sur la paroi n'entrainera qu'une tres faible 
perte de charge, la paroi etant tres mince. 

En effet, 1 'experience montre que la vitesse reelle du fluide est 
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'CR' V CR 



itant in— 



it peine\inferieure a la valour calculie. 

A cause das frottements, la Vitesse r^elle est V r 

CR 

firieur a V_. Le rapport — v — ■ C est le coefficient de vitesse 

C C 

(eompris habituellement entre 0,95 et 0,99). Ainsi le dibit volu- 
■ique est igal a i 



S e V CR ' C « * V CR" 



C C I 

v c 



C itant le coefficient de dibit de l'orifice i 
c 



(3) 



G 




▼ q 

4) Application i me sure du dibit d'une poape i 

Une poape dont le dibit est constant 
aliaente un ricipient ouvert. Le di- 
bit de l'orifice croit au fur et & 
He sure que le niveau dans le rici- 
pient aonte. Lorsque ce dibit est 
igal au dibit de la poape, la niveau 
reste constant dans le rioipient. Ain- 
si la mesure de z M -2 » h permet de assuror le dibit t 



V 



\ 



v/igT? 



itant diternini par un italonage. 



5) Gas d'un gas i Le gas est a la pressioa p t k l'intirieur du ri- 

servoir et se ditend a la pression p # par un petit 
orifice. Si p i n'est pas trop grand de- 







"7* 4 









vant p e , le jet gaseux a un conporteaent 
analogue au jet liquide (contraction avec 
C c - 0,6) aais les effete de pesanteur 
sont nigligeables et le jet est reetili- 



gne. 



Si 



*4fCL, le gas peut Itre con- 



sider conme incompressible. On peut 

de BERNOULLI 8 V? 



done appliquer la relation 

„2 



Pi + 



soit 



et q v -c q 



s V 
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EXERCICES s I) Temps de vidage d'un recipient et loi de variation 

de h avec t. Nous supposons que le rap- 
port S/s est grand pour l'^coulement 
soit conside>e comme p ermanent k tout 
instant. On a V n (t) * ^ 2 g h(t7 
soit q y (t) = C q s ^2 g h(t) ' 
Pendant l'intervalle de temps dt, il 
s'ecoule un vo lume de liqu ide 
dt « - S dh = C q s\j2 g h(t) dt 
soit s dh 

dt - - o q « ymwy 

La durde tjdu vidage est egale a s (on considere que la section 
contracted et 1' orifice ont pratiquement la mime cote) 




P 1 f> s 

h '-)" ■ If* 



dh 
Us^'igh 1 

V 

pour calculer t, il faut considerer la variation de S avec h. 
Si le recipient est cylindrique S = Cte 



*I " 



q 



Sh„ * volume total de liquide 

q 



C s v|2gh 6 '= debit volume k 1' instant initial 
La loi de variation de h avec t est la suivante s 

C. «5-^g*dt 

It 



soit en integrant f h dh . « m . 

J - ?w» = VH* J 0 « 



(\IhJ-\lh 1 ) = C q t 



s 



h = 0 pour t = tj. On retrouve tj = -|- 
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On psut done ecrire i {h*= -"r- ^g 1 ( tj - t) 

soit h = A 2 (t I -t) 2 
La hauteur du liquide est une fonction parabolique du temps. 

II) Form* du jet ; parabola d'axe vertical : r = m)? = m g* 



dans le oas d'un orifice perc< sur la paroi verticals. 




45- = v 



x= V c t 



dz 




L'equation de la trajectoire s'obtient en eliminant t entre les 
expressions de x et z t z g ^2 

c 

Si on considers des jets provenant d'orifice places a differents 
nivsaux, on obtient une serie de paraboles. 

Remarque : Le temps de vidage et la 
forme de la narabole ne dependent pas 
de la density du liquide. Ainsi deux 
volumes d'eau et de mercure s'ecou- 
lent dans le mfime temps. 

Caloul du coefficient de contraction d'un orifice - Orifice de .30RDA 
^ajutage oylindrique interieur.) 

Le liquide est pratiquement au repos tout le 1 onir des parois du r6- 

La lontru. ur L du tube est a6sez 
grand e pour qu'on puisse negliKer 
la Vitesse du liquide au voisina^e 
des parois du reservoir r-t ar.ser 
eoarte^iour que le jet soit libre 
dans le tube (ne recolle pas contre 
les rerois) . 

La repartition des preseions peut done 6tre consideree comr.ie hydro- 
statique. Les forces de pression sur les rarois int^rieures du va- 
se s'^quilibrent deux h. deux sauf en ce qui concerne In section s 
du tube et la section s* projection de s sur la naroi opposee. Les 
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forces de pression exere^es par le liquide sur s* sont egales k 
(ggales aussi aux forces de pression exerc4es par la surface s' de 
la paroi sur 1® liquide) : (p a + ^gh).s 

Appliquons le theoreme d'EULER a la surface S fermee suivante : pa- 
rois du reservoir, surface libre, section droite du jet s^ et sur- 
face laterale du jet jusqu'a cette section droite. 
debit total de la quantite de mouvement sortant par S 

dq T = 2Lf forces appliquees au fluide contenu 
dans S 

CalculonsXF*" : les forces de surface comprennent isui^ggt l'axe Ox 

- forces de pression exerc^es par les parois mouillees 
soit (p a + f gh) s 

- forces de pression qu'exerce l'atmosphere sur la 
portion de jet libre compris entre la surface s 
et la section s du jet : ces forces sont egales 
aux forces de p?ession qui s'exercent sur la pro- 
jection de la surface perpendiculairement a cet 
axe soit i - p & s (voir chapitre pBTATIQUE) 

La somme des projections sur l'axe du tube des forces de pression 

exerc^es sur le liquide contenu a l'interieur de S est done s 

(p a + tgh) s - p a s = £ghs 

En regime permanent (h = Cte), la projection sur l'axe du tube de 

la quantite de mouvement qui sort par S par unite" de temps est 6ga- 

D'apres le theoreme d'EULER : - _ . - . 

C s c 2 gh = C ghs 
On en deduit que le coefficient de contraction C = s /s 

C c = s c /s = 1/2 

Remarque I s En fait pour le debit q il faudrait prendre C sV 

v q c 

C a S /S 

c c' q = C sV = C s V 

H v q c v c c 

et dans ces conditions on a 

C V C C = 1/2 

C q " C v C c = coefficient de debit est egal a 0,§t. C y etant compris 
entre 0,97 et 0,99, C Q est legerement superieur a 1/2. 
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Remarque 2 s Cas d'un orifice en paroi mince 



-n. 



f 



'//// 1 //// 



DU© 



7. 



<J)1T + 7 = ^sgh car R est egal au poids du liquide non hachure" 

(statique) 

ici la composante verticale de R + P est nulle 

Le coefficient de contraction est superieur k 0,5. 

En effet, la vitesse du fluide le long des parois au voisinage de 

1* orifice (portions h.'Churees A et B) ne peut plus Stre negligee 

comme dans le cas de l'orifice de BORDA. 

II y a une zone de depression : la pression est plus foible d'apres 
le th^oreme de BERNOULLI p a gz M = P A + (gz A + f -rfc- 

soit p A = p a + £ gh - £-|^" 

II y a diminution de p A quand est different de zero. 

Ainsi la somme des projections des forces de pression qui s'exer- 

cent sur le liquide sont done superieures a fghs. La projection 

de la quantity de raouvement restant 6gale a f s 2 gh, on a l'ine- 

v c 

gallte^ suivante (th<§oreme d'EULER) 

projection sur axe du tube*-* = £s c 2 gh>^ghs 
soit C c >0,5 
Le coefficient de contraction est superieur a 0,5. 

III-I 3) Ecoulement par dessus un deversoir a bord mince : 

Un deVersoir est un barrage construit en travers d'un ecoulement 
(un canal par exemple) qui perraet de relever le niveau amont et de 
assurer le debit. L» ecoulement peut fetre conside>e comme 

6tant a deux dimensions. Nous pouvons 
appliquer le the'oreae de BERNOULLI. 
Ainsi, la vitesse V d'un point M est 
egale a : 




7-7—7-7 



/ 7 y / 
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Assez loin du deversoir, les trajectoires sont rectilignes 
consequent (cf § 111-1 2) : 



par 



P M , + ^gH = Gte = p a + £gH 0 

D'autre part, la vitesse peut fetre conside"ree comme constante (ega- 
le a V 0 ) le long d'une section normale au sens d' Ecoulement. 
En consequence, il vient : 



M 



+ P 



M 



£gh 



Si le point M est pris a la surface : p M = p & 
2 

-^L = + g(H 0 - h) 

Si on neglige V 0 devant V M , V M = \J 2 g(H 0 -h)' (5) 

La connaissance de V M se d^duit a partir de la mesure de H 0 -h qui 

peut etre effectu^e a. I'aide d'un flotteur. 

Le debit en volume du deversoir est donne par la relation: 

q y = S N[2g(H 0 -h)' 

= coefficient ampirique de dEbit 
S = section dans le plan du deversoir 

EXERCICE Ecoulement par siphon 

Un siphon est une conduite de section constante s qui permet 1 1 e- 
coulement d'un liquide a. partir d'un reservoir, une partie de cette 
conduite se trouvant au dessus de la surface libre du reservoir 

(l ' ecoulement ne s' amorce que si on 
remplit au prealable la canalisa- 
tion) . La section s pourra fetre con- 
sidEree comme suf f isamment reduite 
pour qu'on puisse n^gliger les va- 
riations de cotes dans les sections 
droites ; c'est k dire p Cte. 

La vitesse d' ecoulement V est la 

m 

meme dans toute section droite 




© [G.FANTOZZI], [1975], INSA de Lyon, tous droits reserves. 



- 67 - 



(mouvement eonservatif ) . Le jet a la sortie est sui pos£ cylindri- 
que : ain6i la pression dans le jet est umforme et e*rale \ La 
Vitesse a la surface libre est ndgligeable. 

Appliquons le thdoreme de BERNOULLI entre la surface libre et les 

sections passant par M et- S t 
,2 



t v m 2 
+ e gZ A="V L+ P M + C« Z M =e V m +p a + C^ 



Soit 



V M = 2 g(z A" Z S ) (6) 



^v>thao - S V 



m " S [ 2 «<V Z S>] 



1/2 



p m - P a -^«(v z s> 



La cote du point H n'intervient pas. V croit avec la difference 
z^-Zg mais il existe une limite qu'il n'est pas possible de depas- 
ser. 

En effet, la pression est minimale au point H : 

La limite inferieure de p H est egale a la pression de vapeur natu- 
rente du liquids i p yS » Lorsque par accroissement de z H ~z s » P,, at- 
taint P y g» il forme en H des bulles de vapeur : c'est le phenomene 
de oavitation . Les equations des fluides homogenes ne s'appliquent 
plus et le siphon peut se d£samorcer. 
La difference de cote limite est done : 

^H-^lim - <Pa ~ Pys>/C« 



III-I 4) Pression d'arrSt. Tube de PITOT : 

P ression d'arrfit > Considerons un obstacle immobile dans un fluide 
incompressible en mouvement permanent. 

L'obstacle ecarte les lifrnes de 
courant et en amont de l'obsta- 
//' / cle il y a une zone reduite ou 

Ot>^h\cle la v it esse d'ecoulement est fai- 

ble (ex t pile de pont). 
II existe done une lisrne MR qui 

V\ v \ U \ \\ 
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s'arrete sur 1' obstacle en R. Rest un point d'arr&t ; V R 0. 

Appliquons la relation de BERNOULLI : 

,,2 

P R + C« Z R = P M + f gZ + C 

V 2 p» p * v 2 

P R = P M + C-2— 0U Tf— = Ti™ + ~2T 

Si on place un tube pi^zometrique en R, la cote du fluide donnera 
la quantity ~^g~' 

Si en M on met une prise de pression statique, le tube indiquera 
if 2 

— . La d4nivellation entre les 2 tubes h est 6gale a . 

Ze r~ t=t- ' 2g 



On a par consequent 



Sans obstacle, la ligne de courant MR se poursuivrait et en S' la 

vitesse V„, serai t different© de zero . 
rt 

D'apres la relation de BERNOULLI : 

f V M ^ V R< 

*-g=- + P M + £»z = + P R ,+^gz R 

Si le point M n 1 est pas trop voisin de R af in qu'on puisse admettre 
que l'ecoulement n' est pas modi fie *par 1 ' obstacle (en M X ! : la quan- 
tity pr^c^dente est aussi egale a p R + f> gz R en presence de 1 1 obsta- 
cle. On en d^duit done : 

p R(+ ^—j- = p R 

La pression au point d ' arret ou pression d ' arret p R est tou jours 
supe>ieure p R , i c 1 est la pression totale de l'ecoulement preexis- 
tant en R' (pression statique +pression dynaraique) . La prise de 
pression R s'appelle pour cette raison prise de pression totale . 

REMARQUE t Gas des gaz : On peut negliger ^gz devant p 
On a ainsi p R = p^ + £ — 



V .. P ^" Pm ' ( 9 ) 

Les pressions p R et p M peuvent Stre mesurees \ 1 ' aide de raanometres 
ou p R -p M a l'aide d 'un manometre dif f£rentiel . 
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Tube de PITOT simple t C'est une ante me de pression totals. Ses 

dimensions transversales doivent fttre assez 
faibles afin de limiter les perturbations. Le tube recourbe' relie* 
au manometre remplace 1' obstacle. 

double $(ou tube de PITOT-PRANDTL) 
Appareil qui permet de me surer las vi- 

-s ~~ ,,f tv>ano- 

{< i „, e vv-t tesses d'ecoulement en un point. 

II est constitu£ d'une prise de pression totals R et d'une prise 

de pression statique M relives a un manometre differential. 




L'appareil doit Itre placl parallelement aux lignes ds courant. 
La vitessa d'scoulsstSnt est donn^ par les relations (8) ou (9). 
Cet appareil est tres utile oar il permet un* me sure facile des vi- 
tesses 6ooulement (ex t vent, fluide, etc...). 



III-I 5) Phenomene de VENTURI - Mesure des debits t 
Phenomene de VENTURI i 

Considlrons une conduite dont la section S est variable. 

Si on place 3 prises ds pression 
statique a l'entree, a l'sndroit 
Is plus etroit st a la sortie, on 
observe des niveaux differents 
dans les tubes piisoaetrlques. 
Si on admet que la conduite se com 
porte comme un filet de courant (1 
vitesse V est uniforms dans une 
section quelconque), on psut appli 
quer le th6oreme ds BERNOULLI t 

ear * ^ " ° te 
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D' autre part, la conservation du debit volumique s'lcrit t 



q v - Sj Vj - S 2 V 2 - S 3 V 3 = Cte 



Si S diminue, V augmente et par suite p diminue. 

JUL 

S et p variant done dans le mftme sens. Ainsi dans la parti® la 
plus Itroite de la conduite, il existe des depressions qui peuvent 
fetre importantes. 

APPLICATIONS t I) Trompe a eau (vide) t Supposons qu'un point A, la 



2) Carburateur d'un moteur a explosion t L 1 aspiration et la pulve- 
risation de 1' essence sont provoqules par 1 ' Icouleraent <Je I'air de 
carburation issu de 1 ' atmosphere dans un conduit qui coraporte un 
r<H x ecisseraent ("bus*") au niveau duquel 1 ' accroissement d® Vitesse 
entraine une depression qui aspire 1' essence. 

3) Mesure de dibit t tube de VBNTURI j 

C'est une canalisation convergente-divergente qu'on interoale dans 
la conduite dont on veut mesurer le dibit. 




/ 



pression absolue de 1 'eau est la 
pression atmospherique. 
Au niveau da l'ltranglement, la pres- 
sion p B est inflrieure a p A , AJnsi 
l'air se trouve aspirl en B et en- 
tralnl par l'eau. 
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Le d^bit de la conduite est egale k 



1 



T-Po 



Pj-P 2 



2 pl 



(10 bis) 



(10) 



se determine a l'aide de deux tubes piezometriques 
ou a l'aide d'un manometre dif ferentiel . 
Si le tube est dispose horizontalement ou si le fluide est un gaz, 
la me sure de p j - p 2 se ramene & la me sure de pj- p g . 

REMARQUE « I) Le divergent ne doit pas §tre trop ouvert pour que 
lea lignes de courant ne se decollent pas de la paroi (<I0°). 
Le convergent peut avoir un angle plus grand. 
2) L« d<bit reel q yr = £q y 

"Jetant un coefficient correctif qu'on determine par un ^talonnage. 

EXERCICE ! Mesure des debits i diaphragme (organes d^primogenes plus 
simples que les tubes de VENXURI mais moins precis). 
Le diaphragme est une ouverture circulaire effectuee dans une paroi 
mince. La section de la veine fluide pr^sente un mini-urn lege>ement 



— - . 

JL 



en aval du diaphracrme. Calculer 
le debit de la conduite en fonc- 
tion de h, d et D diamotres du 
diaphragme et de la conduite. 
Appliquons le th^oreme de BERNOUL- 
LI i 



e> c 



A 

2p 



2ir 



. PC 



S A V A = S B V B 



(les cotes sont identiaues) 
= S C V C 



D' autre part 

La valeur de est legerement inferieure a Sg . Sunposons en Iere 
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approximation que S B = S c , V B = Vq. 



h n kj \j 



Soit 



p A- p C 



h = 




En fait cette relation n 1 est VrTgoureuse car la section reduite n'est 
pas £gale a S c et parce qu'il y a des pertes de charge. 
On introduit par consequent un terme correctif : 

Q vr£el = ^ Q vtheorique ^~ 0,5 

Tuyere : Instrument de mesure intermediaire entre le diaphragme et 
le VENTURI. 



Le fluide epouse le contour de la 
tuyere mais il n'y a pas de partie 
divergente. La relation donnant le 
debit est la meme que pour le VEN- 
TURI ou le diaphragme : 



v B 



2gh I 



I- 



tll-I 6) Representation graphique de l'ecoulement 
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On porte sar la verticale passant par le centre de gravite" d'une 
section : 

- une longueur OG^. = Zj = cote du point Sj qui represente le terme 
de pesanteur. 



cote du point Pj est la hauteur piezometrique (niveau du tube pie- 
zometrique relie a une prise de pression statique dans la section 
Sj). Quand on considere les diff ^rentes sections de la conduite, 
le point P T decrit une-courbe qu'on appelle ligne piezometrique . 



tique. 

Le point Qj d6crit la ligne de charge . La cote de Qj est 6gale a. la 
charge totale dans la section Sj. 

Dans le cas des fluides parfaits, s'il n'y a pas de g£nerateurs ou 
recepteurs, la charge totale est constante et la ligne de charge est 
horizontale (plan de charge). Pour ies iluides r^els, il y a tou- 
jours des pertes d'^nergie de sorte que la ligne de charge est des- 
cendante. La difference represente la perte de charge. 

II-2 Applications du theoreme d'EULER : 

II-2 I) Action exerc^e par un .jet sur une plaque plane fixe s 

Considerons une plaque plane faisant un angle^avec un jet horizon- 
tal de Vitesse uniforme V. Le fluide est devie par 1 ' intermediaire 



Appliquons le theoreme d'EULER a la surface ferm^e S delimitee par 

les 2 sections I et 2 et la surface exterieure du jet. 

Les forces exterieures qui agissent sur le fluide contenu a. l'inte- 



- une longueur ii T P T = 



qui represente le terme de pression. La 





de 1'epaulementE. 

Les vitesses dans les sections I et 
2 sont identiques 5 en efiet l'ecou- 
lement etant a pression constante, 
1 ' application du theoreme de BER- 
NOULLI conduit a l'^galite suivante : 
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rieur de la surface considered sont : 

- le poids du fluide i vertical (on neglige les effets de pesanteur) 

- les forces de pression : p a pour les elements en contact avec 

1 11 atmosphere et p +p pour les elements en contact avec la plaque. 
La resultante des forces de pression atmosphferique est egale a la 
poussee de l'air et par consequent « 0, 
^Pour les sections droites I et 2, on a : 

p + { g z = Cte = (pf (g Z ) bQrd jet = p a ♦ e«Wd jet 
Si on neglige les effets de pesanteur, on a p = p a . j 
Seule intervient la resultante des forces de pression p. 
La resultante des forces de pression de la plaque sur le fluide est 
done £gale a (theoreme d'EULBR) : 

R ss q V (n 2 ~nj) n 2 , n-j £tant les vecteurs unitaires 

perpendiculaires aux sections I et 2 
= V n^ V 2 m V "n 2 

L 'action du jet sur la plaque est <5gale h $ 

if = -T= q V (nj-nj) (II) 

La composante de cette action a pour valeur : 
— * — * 

Q x = Q.nj = q V (I-coso( ) 
(n^."n^) = cos<L 

Q x = q V (I-cos«0| (12) 



Ce r^sultat est ind^pendant de la forme prise par l'ecoulement entre 
les 2 sections I et 2. Le d^bit massique q etant a pVs (s etant la 
section^ droite du jet). On a : ' 

Q x = £ V 2 s (I-cosV ) 

La poussee suivant la direction initiale du jet se fait dans le sens 
de l'ecoulement. Elle est maximale pouro^=1Pet vaut alors 2{£V. 
Remarque : Ce phenomene est utilise dans les turbines hydrauliques 

PELTON mais dans ce cas l'auget de turbine n'est pas fi- 
xe. Si V 0 est la vitesse de translation reetiligne uniforme paral- 
lele a u jet initial. 
© [G.FANTOZZI], [1975], INSA de Lyon, tous droits reserves. 



- 75 - 



i ec a 

l'auget, le m&me raisonnement que precedemment conduit a : 

Q x = V q" (I-costf ) = £ V ,2 s (I-cos^ ) 
V etant vitesse relative du jet, la vitesse re"elle etant egale a 

v - v + v 0 

soit Q x = £(V-V 0 ) 2 s (I-cos«C) 

Exemple % q=3o l/s diametre du jet 3 cm c(= 165° 

Q x = 2500 N = 250 kgf 

III-2 2) Action d'un fluide sur une tuyauterie : 

Conside>ons deux sections I et 2 d'une conduite reliee a un r£ser~ 

2-">&f servoir dont la surface li- 

I z =o 



^ _ 0 bre est prise comme plan 

de reference. 



, ft* 



Appliquons le theoreme d'EULER : 

q vt - q vT = *R h + 1p*^ , . 

2 I ^ forces volume (pesanteur) 

forces surface 

qV*£tant applique au centre de gravitfi de la section 
R comprend : 

- les forces de pression sur les sections I et 2 exercees par le 
fluide soit f ,„ -* , f . jo ^ 

ysj * js. 



J 2 

- la reaction R^, des parois du tuyau sur le fluide (reaction egale 
et opposee a la force Q exercee par le fluide sur les parois du 
tuyau) . 
On a done : 
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q ( ^2"^I ) = J P I dS I "l " J ?2 dS 2"^> + P*+ % 

d'ou 0*= -"r£ * q (V^-V^) + Jpj dSj -jp 2 dS g n^ + P* 

Nous allons calculer successivement 1' action statique du liquids 
sur les parois puis 1' action due au mouvement du fluids. 

Cas des forces de pression sur les sections I et 2 i 

Pour un point de cote z, la pression relative dans un fluide au re- 

pos est p 1 = - £ g z (voir figure). 

Soit p la pression au mfeme point dans le fluide en mouvement. 

On a : JH 0 nrr „ , 

P = p + rgz = p-p' 

soit p = p + p' 

Ainsi la poussee f p T dS T est la somme t 
K 

- d'un vecteur repr^sentant la poussee hydrostatique sur Sj (terme 
p') et applique i- au centre de poussee Q T 

- d'un vecteur (terme p**) t 

Dans la section droite Sj, p j = Cte done Pj est applique* au centre 
de gravity de Sj. 

P* - P * Sj nj 
11 en est de mfeme pour la section 2. 

Cas de la reaction R T J 

peut fttre decompose* en deux parties : R^ et R*^ 
— * — ■» i 

Q' = - Rj, = poussee hydrostatique du fluide sur les parois du tuyau 
Qt^ - - R^* = poussee suppl4mentaire due au mouvement 

On aura ainsi t _± -» — » -» » -si* » — s» 

-(Rj + R*) = q (Vj-V 2 ) + P+ P£+P*+P'+P» 

Dans un fluide au repos, I 1 ensemble des pouss£es qui s'exercent 

sur une surface fermee S est eg ale et opposes au poids du fluide 

contenu a l'interieur de cette surface s 
Ainsi P* + PJ + + "r£ = 0 Q' = T+ p|"+ P^ 

d • ou - RjjT- q ( - Vg*) + P^ 
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La pouss^e supplementaire due au mouvement du fluide s'exercant 
sur la tuyauterie est done egale a : 



Q*= (p* Sj + q Vj) nj. - (p* S 2 * q Vg) n g 



(13) 



La poussee hydrostatique augmente avec la cote z alors que la pous- 
see dynamique n'en depend pas : elle depend de la vitesse du fluide. 

Exemple s Considerons une conduite coud^e a 90° de diametre int£- 

rieur 100mm, avec un £coulement de 30%/s, la pression ef~ 

fective (ou relative : p^ = p-p 0 ) du fluide £tant de 20N/cm „ 

e a 

Calculer la r^sultante des forces dues au mouvement du fluide ( on 
neglige le poids du liquide ) 

S I =S 2 = 78,5 cm2 




p i s i =p 2 s 2 

V 



78,5 x 20 = 1750 N 



I~"2 



V„- ™ I" \ K =3 ,82 m/B 



78,5 10 



qVj=qV 2 = 30 x 3,82 = 114,6 N 



1 



= 1000 kg/ 3 



Soit (1750 + 114,6) ^fiT 1 

*f 2360 N 

I litre = IO 3 cm 3 = IO~ 5 m 3 



III-2 3) Reaction d'un .jet issu d'un reservoir s 

Consid^rons un reservoir perc4 d'un orifice. Nous connaissons la 
repartition des vitesses dans le jet et le reservoir (application 
du theoreme de BERNOULLI). 

Appliquons le theoreme d * EULER a la 
surface ferm^e S indiqu^e sur la fi- 
gure (On conside>era la projection 
horizontale) . 
On a *' = q V 

F = resultante suivant 1' horizontale des forces applique'es au fluide 
contenu dans S (forces de pression exerc^es par les parois du re- 




servoir) 



r = (s c v V = (s c V2 



(14) 
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Inversement, le reservoir subit une poussee Q £gale et opposee qui 

est la reaction du .jet . n 0 „2 
" — Q a ~ \ B c v 

Si le fluide est un liquide, on peut appliquer le th^oreme de 

TORHICELLI i Q a ^ ^ y 2 m g ^ ^ g z (J5) 

la reaction du jet ne depend que du liquide contenu dans le reser- 
voir. Cette reaction tend a deplacer le reservoir en sens inverse 
du jet, 

III-2 4) ££MS£S d'un reacteur ou d'une fusge : 
a) Cas d'un reacteur i 

La propulsion est obtenue en ejectant a l'arriere les gaz produits 
par combustion du petrole dans de l'air prealablement comprime. 
Cans le cas d'un mouvement de translation uniforme de l'avion, l'air 
•ntre par le reacteur a la vitesse relative V par rapport au reac- 
teur (sysfceme de reference lie' a l'avion), les gaz etant ejectes k 
la vitesst relative W. 

On peut neglige r le debit massique de petrole par rapport a celui 
de l'air | le debit massique des gaz Ejectes q est done egal k ce- 
lui de l'air entrant. _ 

Appliquons le theoreme d'EULER *~ ~ ~ "^-c^:-, 

l — > -i ^ 

k la surface fermee trac4e en I \j — 1 ^j* 

pointing: ~ - ' " 

La re suit ante des forces de pression \ 

est nulle ? la poussee sur le reacteur est egale a : 

Q = q (W - V) (16) 

Cette poussee depend de la vitesse V de l'avion. 

Si V = 0, Q = qW : la poussee est maximum. 

La vitesse maximum theorique de l'avion serait la vitesse W d 'ejec- 
tion des gaz ; dans ce cas Q = 0. 

En mouvement uniforme, la force Q = q (W - V) est compensee par la 
force de frottement de l'air s'exergant sur l'avion. On peut calcu- 
ler le rendement de la propulsion. 

La puissance employee k cette propulsion est egale a : 
Q V = q V (W-V) 
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La puissance fournie a l'air pour accroltre son energie cinetique 
de 1/2 q V 2 a 1/2 q W 2 est egale k : 

P = 1/2 q W 2 - 1/2 q V 2 = 1/2 q (W+ V) (W-V) 
D'ou le rendement de la propulsion s 

R _ QV 2V_ _ 2 , T7 . 

R " -p W+V~ " TTTT (I7) 

v~ 

Le rendement est toujours infe>ieur a I (=1 si V tend vers W). 
b) Cas d'une fusee : 

Une fuse'e contient a la fois le combustible et le carburant. 

La poussee est egale a : 



-c 



y [ ^ ^."V^ Q = qW (18) 

f~ Cette poussee est independante 

de la vitesse V de la fusee. 



La puissance mise en jeu est : 
QV = q WV 

Dans le vide, la fusee subit une acceleration & egale k : 

M = masse de la fusee a 1' instant t 
Cette masse diminue la quantity q par seconde si bien que, si ^ est 
constant, 1 'acceleration augmente sans cesse jusqu'a ce qu'il n'y 
ait plus de c«mburant. 
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CHAPITRE IV 

IV EQUATIONS DYNAMIQUES DES i-'LUIDES VISQUEUX 



IV- I Coefficient de viseosite : 

IV-I I) Definition : 

La notion de fluide parfait n'est qu'une premiere approximation 
quand on considere la dynamique des fluides reels. 
Alors que dans les fluides narfaits, la tension qui s'exerce sur 
un element de surface dS est tou;jours normale a dS, dans les flui- 
des reels les actions de contact presentent une composante tangen - 
tielle . qui joue le r&le d'une force de frottement. Ainsi, une 
nouvelle propi£t£ des fluides est introduite : c'est la viseosite 
ou frottement interne . 

La tension s'exercant en un point M 
peut se decomposer en deux parties : 

- composante normale : pression 

- conroosante tansrentielle 
La viseosite introduit la composante 
tangentielle et modifie la composante 

normale. Ces forces tangentielles dites forces de viseosite s'op- 
posent au glissement relatif des couches voisines. 
L ' experience de principe suivante permet de montrer 1" existence de 
ces forces de viseosite' : 

Soit une couche fluide assez mince comprise entre plaques planes 
paralleles ; l'une A 0 B 0 est fixe, l'autre AB se deplace dans son 
plan avec une vitesse V 0 . 

On constate que le fluide en contact 
avec la plaque A 0 B 0 est immobile alors 
que le fluide en contact avec la plaque 
AB a la vitesse V 0 (adherence aux 
plaques) . 

La vitesse V du fluide en un point M passe progressivement de 0 a V 0 
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quand z varie de 0 a h. 

Si h et V 0 sont assez faibles, les couches fluides paralleles aux 
plaques glisaent les unes sur les autres sans se melanger (le mou- 
veraent est dit laminaire et noue preciserons ulterieurement les 
conditions pour obtenir un tel £coulement). 
On admet alors comme prineipe la relation suivante : 

- - Cte (I) 

soit u = Ctexz ou du = Cte X dz 

On a done, suivant la direction z, un gradient de vitesse constant 

-SS- - -fr - •« »> 

Considerons une couche infiniment mince de fluide, d'epaisseur dz. 



La couche superieure exerce sur 
cette couche une contrainte tan- 
gentielle +t?(les vitesses des 2 

Y ~— » ~C couches etant differentes, elles 

^ z A ~£?ZZZ22LZ ^2L /////// subissent des actions de viscosite, 
~t actions qui tendent a reduire les 

differences de vitesse : la eou- 
^ q che superieure tend a entralner la 

couche considered). 
La couche inferieure exerce une contrainte • 'C • Ainsi la couche est 
soumise a une force suivant l'axe Ox nul (prineipe d'inertie : 
F = m?, 2T= 0). 

Ainsi la plaque AB exerce sur le fluide au contact une contrainte 
+TJ alors que le fluide va exercer sur A 0 B 0 une contrainte 
La loi de NEWTON (1713) indique que V est proportionnel au gradient 
de vitesse ; on pose 



du 
ET 



(2) 



^Metant le coefficient de viscisite' dynamique (coefficient de frot- 
tement interne). Ce coefficient depend de la nature du fluide, de 
sa temperature et de sa pression. 
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Remarques t I) Dans le cas de l*ecoulement laminaire, seule existe 



Les unites sont : 

-systeme C.G.S. : poise (P 0 : g cm -1 s -1 ) 
-systeme MKSA s Poiseuille (PI : H s m -2 ou kg m~ I s~ I ) 
I PI = 10 poises 

3) Fluidite : 

Elle est definie par 1/^ 

L'unite" C.Cr.S. de fluidite est le rhe. 

IV-I 2) Fluides NEWTON IKNS i 

Un fluide est dit NEWTONIEN si le coefficient de viscosite )k ne de- 
pend que de la pression et la temperature (ne depend done pas du 
gradfnt de vitesse ^) . 

Si/l depend du gradient de vitesse, le fluide est dit non NEWTONIEN. 
Dans ce cas, le graphe (t» , -j^-) n'est plus une droite. 



la contrainte tangentielle de viscosite. 



Pour I'ecoulement turbulent, au terme 5 |H s'ajoute un terme 
supplementaire . 

2) Unites- : La dimension deWest (relation 2) 




T f 





fluide NEWTONIEN 




Fluides non NEWTONIEN 
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Les liquides pUrs et les gaz sont HEWTONIENS. 

Dans la suite du cours, nous ne considererons que les fluides 

NEWTONIENS. 

IV-I 5) Coefficient de viscosite cin^matique : 

II est defini par^= ~^r~~t f extant la masse volumique du fluide. 
I' Equation aux dimensions donne 
-Iwm-I 



L MI 



= lV 1 



ML 
2-1 



D'ou les unites : 

2 -I 

-systeme C.G.S. : cm s ou stokes (St) 

2 —I 

-systeme M.K.S.A. : m s ou myriastokes : maSt 
4 

Elle vaut 10 stokes. 
Le tableau suivant donne quelques valeurs des coefficients de vis- 
cosity dynamique et cin^matique de quelques fluides k 20°. 



LIQUIDES ^ZO°C) 




GAZ 


(jio'c) 






>(P1) 


Itasi 




^ (PI) 




C (°K) 


Eau 

A i?nflique 
Benzene 

Acetone 

Mercure 

Slycerine 


IO" 3 
1,2 IO -3 
0,65 IO -3 
0,35 IO -3 . 
1,55 IO" 3 
800 IO -3 


' I 
1,5 
0, 74 

0,115 
635 


Air 

co 2 
ci 2 

H 2 
°2 


1,8 IO -5 
1,45 IO -5 
1,35 IO -5 
0,9 IO -5 
1,2 IO -5 


15 
8,1 

109 
15 


114 
277 

70 
138 


Hui rfc$n 


1300 IO -3 












^rre uooo 


10 (ve IIIuide) 












^erre 5000 


IO 13 













R emarque : La viscosite des melanges n'est pas donne'e par la regie 
d'additivite. 

Ex : eau et alcool a 0° C : 64 °/ 0 eau et 36 % alcool :JJ{ = 7,3 I0 -3 P1 
(1,8 I0 _3 )(I,8 IO -5 ) 
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IV-I 4) Influence de ia pression et de la temperature sur le coefficient 
de viscosite : 
I - Gas des liquides s 
a) Bffet de la temperature t 

La viscosity de tous les liquides decrolt quand la temperature augmen- 
te. Aucune loi g^nerale dejA = f(t) n*a pu 6tre donn^e. 




On peut repr^senter la loi de variation par une formule empirique 
de la forme : u 

j{ JL — 

eau = 1,8 I0~ 5 P1 

<^ = 0,0337 f> = 0,00022 
Une theorie de la viscosite due a ANDRADE donne pour^Uune variation 
du type ; A g B/T A et B etant des constantes 

(formule d ' ANDRADE) 
Cette loi est bien verifiee pour un grand norabre de liquides. 

b) Effet de la pression : 

La viscosite augmente avec la press-ion pour tous les liquides mis a 
part l'eau. 

La loi de variation est en general exponentielle : 

A. m |p/po)-3 

Dans le cas de l'eau, les variations sont plus complexes : 
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4000 



gooo 



2 - Gas des gaz 1 
La th^orie cin^tique des gaz permet d'expliquer 1' existence de la 
viscosity. Elle nous donne 1* expression suivante t 

kfcl (3) 

k = 0,499 0" = vitesse moyenne des molecules 

de gaz. 

^ = masse volumique du 1 = libre parcours moyen 

gaz 

Le libre parcours moyen 1 est £gal a (theorie cin^tique des gaz) t 
1 - I 

"nyi'r 2 n 

-<J"6tant le diametre de la sphere d' action mol^culaire (distance k 
laquelle se trouvent les centres de 2 molecules au moment de leur 
collision) . 

- n ^tant le nombre de molecules par unite" de volumes. 
La vitesse c" est egale a ; 

m = masse d'une molecule 
M = masse molaire = N m 



kT 
m 



mn 



pV = RT 
pour I mole V * M/|» 



p/{ - RT/M 



Ainsi le coefficient de viscosity est egal k 1 

(4) 



/- 



0,112 m c 
1 "2 



D'apres cette relation, la viscosite d'un gaz est ind4pendante de 
la pression et vane comme >|7 . 

Pour les gaz reels, on observe que^lvarie peu avec p. Dans le domai- 
ne des fortes pressions (supirieures a quelques dizaines d 1 atmos- 
pheres) , la viscosity augment e avec la pression. 
Contrairement aux liquides, la viscosity des gaz augment e avec la 
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temperature. 

D'apres la relation (4), est proportionnel a vfTl En fait 
SUTHERLAND a etabli une relation plus exacte de^J en tenant compte 
des actions intermoleculaires : l+j£L 

C coefficient caracteristique du gaz : Constante de SUTHERLAND 
(en fait C crott avec la temperature legerement). 
Exemple s air : 20°C / = 1,8 I0~ 5 PI 

IOO°C ^4=2,2 IO~ 5 PI 

300°C Jl = 4,4 IO" 5 PI 
Remarque s : I) la viscosite" des gaz est a 0°C de l'ordre de I0~ 5 P1 
soit 100 fois plus faible que celles des liquides a 0°C (de l'ordre 
de IO" 3 ). 

2) la viscosite cinematique des gaz est en general superieure a cel- 
le des liquides. Elle augmente naturellement avec la temperature 
alors que celle des liquides diminue. 

La viscosity cinematique des gaz est inversement proportionnelle a 
la pression, la masse volumique lui etant proportionnelle. 

IV-2 Equations dynamiques des fluides visqueux : 

IV-2 I) Relations entre les contraintes et les vitesses de deformations i 

Les relations generales entre les contraintes et les vitesses de de- 
formations ont ete donnees par NAVIER en 1823. 

Le tenseur des contraintes d'un fluide visqueux peut 8tre determine 

a partir des considerations suivantes. 

Dans un fluide parfait, le tenseur est = - pSj_j 

p etant la pression donnee par liquation d'etat du fluide,, 

Pour un fluide visqueux, on ecrira z 

<r la = - P % u + # ti (6) 

. etant le tenseur "visqueux" des contraintes. 

Les processus de frottement interne dans- un fluide ne se produisent 
que dans les regions du fluide ou il existe un gradient de Vitesse. 
AinsiO*' ^ doit adpendre des derivees de^vitesse par rapport aux 
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coordonne'es. 'Si les gradients de Vitesse ne sont pas trop grands, 
on pourra admettre que les contraintes <f' . . sont des fonctions li- 



pression detT' i; . puisque les c 
s'annuler si la vitesse iTest 



n^aires des derivees premieres de la vitesse (hypo,these de NAVIER). 

i 

II ne doit pas exister de termes independants de -^-jj — dans l'ex- 

contraintes fie viscosite^"' ^ doivent 
constante. 

LesO i; j doivent aussi s'annuler lorsque le fluide est soumis a. une 
rotation d' ensemble uniforme, le frottement interne dans le fluide 
etant alors nul. SiSlest la vitesse angulaire de rotation uniforme, 
la vitesse v d'une particule est 6gale a. : 



Les sommes 



En effet 




?=0M 



■^1 



s'annulent si v =-S.a r 



+ Txf" 



= 0 



Vi et 



^u. 



^""' . . contiendra done des combinaisons symetriques des termes — i-. 

Remarque : L'hypothese de NAVIER n'est qu'une approximation, veri- 
fiee par l'experience pour des vitesses de deformation 

faibles, 

Dans le cas d'un fluide isotrope, le tenseurtf*'.. satisfaieant aux 
conditions precedentes pourra prendre la forme suivante : 

T57 



"*u. ~ du ~ ~ 
X 1 1 ^ X i ^ X e 13 



(7) 



Jk et X etant des coefficients independants de la vitesse qu'on ap- 
pelle coefficients de NAVIER n' est autre que le coefficient de 
viscosit4 dynamique deja defini etAie deuxieme coefficient de vis- 
cosity. 



^u 



Remarque J La quantite 



W + Ty~ 



^v 



d^r 



est la 



vitesse de dilatation cubique 0 (<T"= volume massique) 
•*u„ 



e = 



ee 
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On peut done ecrire 



<•*..= A6<?. . + 2 li v. . 

1J 13 A 10 



(7 1 ) 



v.. £tant la vitesse de deformation (cf. chap. II) 
l * ^ u. ^u. 

v ij " ^ (TxJ- + 
La contrainte totale^.. est egale a 



du . 



ou 



"<r n = - P+ A6 + 2^v n 
<T 53 = - p 2^v 53 



' (8) 




23 = 2 /* V 23 



La force l^qui agit sur une unite d'aire de surface, "n*etant le vec- 
teur unite normale a cette surface, a pour composantes : 

p i " +<r ij n 5 = ( -P n i + 2 / v i 3 n j 

Le terme -pn.^ represents la pression usuelle du fluide. 



Les autres termes representent la force de viscosite agissant sur 
la surface. 

La composante suivant la normale a la surface ou pression est mo- 
dified : 



P. n. = (-p+A© ) + 2^v. ^ n. n 
IV-2 2-) Equations de HAVIER-STOKES : 



Reprenons les equations dynamiques (voir cours 3GP et chap. I M£- 
canique des fluides) 



Remplagons ^ . par l'expression (8) 



On a : 



^ tX ft _ \ \ © 1 



ij "4x 

J 



(Ae.p) =A^_ - 



ax: 



^ X? «x . &x 

o i ; 
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}> u. Su. 
2 v. . = r-=— + v 



Les coefficients A eyl etant consideres comma constants (dans la 
plupart des cas, leur variation dans le fluide est tres faible). 



Or on a : 



Vu. 



A- =Au. 



Vu, 



9= div V = e 



^x e 



On obtient ainsi les Equations de NAVIER et POISSON s 

du. 
i 

-ar— 



"P T#- = X i + ~f~ Ti~ + ^" Au i" ~dT" 
-i- gra3*p = X*+ grad divT+^AlT — 

c " dt 



(9) 
(9') 



Dans le cas d'un fluide incompressible , div V=0. 

Liquation du mouvement d'un fluide visqueux incompressible s'ecrit 
(equations de NAVIER- bTOKES) : 



du. 

ir- 



on peut expliciter 



du. 

"cfF" 



]5 

+ U . <<r i 

3 d x 



(10) 



( 



TT 



0 pour un mouvement permanent) 



i 



sous forme vectorielle : 



- -B- graci p + "f*+ $AV* = 

V J \ ^ 

forces de pression vol 



dV 
TT 



(10') 



Slume 



viscosite 



Remarque t Si^JJ = 0, les Equations de NAVIER-STOKES se reduisent aux 
equations d'EULER. 



V 
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Equations intrinseQaes : 

On peut rapporter les equations dynamiques du mouvement a des axes 
lies a la trajectoire d'une particule. 

axe Ox est la tangente"T a la trajectoire 
A M de la partieule M. 

axe Oy est la norraale principale n dirig^e 
vers le centre de courbure C. 
On considerera un £coulement permanent 
done ^ = 0 




Dans le cas d'un fluide incompressible dans le champ de pesanteur, 
on a : f = Cte et F = - gracf gh 

Les equations de NAVIER-STOKES s'ecrivent alors : 

= - ^grad gh - graS p +^&V* = - grad (p+£ gh) +)<Av* 

— ^ 

Projetons cette egalit4 vectorielle suivant la tangente t et la nor- 

-* dv 

male n, Les composantes de ^ sont respectivement : 



dV 4 



= 

n 



dt 



TT 

Tt~ + Ts" 



ds 



dS _ " n 
~dl VF 



+ V 



t TT 



t TF 




On a R dt( = ds 



dc^ 



R = 



<1V_ 



ds ds 



~3* — W r 



Pour un ecoulement permanent, -j-j— = 0, d'ou : 



(ID 



dV 



(V^ et V sont identiques) 



if = 

n 



R 
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Ainsi pour un ecoulement permanent, 1« equation de NAVIKR-STOKES 



s'Scsrit 



-^-(p+f «h) + (^tify 



i 

Multiplions les 2membres de la premiere equation par ds = V dt 

(V dV = - d(p+fgh) + (JK^v) t ds 

en restant sur la trajectoire 

so it. v 2 

pd(JL_) = - d(p+^gh) - 

2 

dtfe'tant le travail fourni suivant la trajectoire par les forces de 
frottement (par unite" de volume de fluide). 
On peut encore 6crire : 



K -ftg - + p + ^ crh) = - dtT 



(12) 



d"5Test un travail m^canique perdu par le fluide (travail transfor- 
ms en chaleur par frottement visqueux). 

Si d 1> ■ 0 (^A = 0), on retrouve 1'dquation de BERNOULLI du fluide 



parfait incompressible 



+ p + ^ gh = 



= Cte 



le long de la trajectoire. 
Remarque : cas des trajectoires rectilignes paralleles : Dans ce 
cas, la composante j}*&V) n est nulle'. Comme R est°o, 

«^j-(p + £gh) = 0 soit p + £ gh = Cte (15) 

Dans un plan perpendiculaire aux trajectoires, la repartition des 
pressions est hydrostatique, comme pour un fluide parfait. 
Aux equations de NAVIER-STOKES, il faut ajouter l 1 equation de con- 
tinuity (equation H->), 1' Equation d'Etat du fluide et le princi- 
pe de la conservation de l'energie, soit au total 6 equations pour 
determiner 6 inconnues u, v, w, p,^, T. 

Condition aux limites s On admet d'apres 1' experience que les ele- 
ments fluides au contact d'un solide adherent au solide sans glisser . 
Le long de la paroi, la vitesse du fluide est identique a celle de 
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la parol. De ce fait, la limite du mouvement d'un fluide visqueux 
dont le coefficient de viscosite tend vers 0 n'est pas en general 
un mouvement sans discontinuite de fluide parfait. 
Remarque : Les Equations de NAVIER-STOKES ne sont pas lineaires et 
de ce fait n'ont pu Stre int£gr£es que dans des cas partieuliers. 

IV-2 3) Dissipation d'4nergie dans un fluide incompressible : 

La viscosity entralne une dissipation d'^nergie qui est transfor- 
med en chaleur. Considerons le cas d'un fluide incompressible. 
Le travail de deformation (reel) par unite" de masse est egale a : 

T Vn dt = "T V^d = ' i~ Viidt + i~ < ^ dt 
= - ^-6dt + j-*"^ v 4j dt (14) 

D'apres le premier principe de thermodynamique (principe de con- 
servation de l'^nergie), on a la relation suivante : 

dW + dQ = dV + dC (15) 
V etant l'lnergie interne du systeme 
C Anergic cinetique 

dW travail recu par le systeme (travail des forces ext^rieures) 
dQ quantity de chaleur regue par le systeme 

La chaleur regue par le systeme dQ est egale a (second principe) 

dQ = T dS - df (16) 
avec df>0 oudQ>TdS 

df = 0 pour une transformation reversible : dQ = TdS 
df )>■ 0 correspond aux transformations irreversibles (irreversibili- 
ty mecaniques, physiques ou chimiques). Les irreversibilites rae- 
caniques interieures sont les frottements interieurs ou viscosite. 
Le travail reel des forces d'inertie est egal a -dC, le travail t de 
deformation est done egal a dW - dC = dU - dQ 

D'apres le 2eme principe, dW = dW - dC = dU - T dS + df 

def 



dW . . = dU - T dS + df 
def 



(17) 
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Pour un fluide a l'6quilibre, dU - TdS = - p deT(transf. reversible) 

Le Ier terme de la relation (14) qui est £gal a. p d$" Atr - 

~ ■* qr~= - P da 

(tf*= volumique massique) correspond done a des vitesses de deforma- 
tion infiniment petites. Dans ce cas, la transformation est rever- 
sible et df = 0, Le Ier terme repr^sente done dU - T dS alors que 
le second terme represente le travail dissip4 df . On a done : 



fdf jVi3 dt ■ 0dt 



(18) 



fiest appelee fonction de dissipation. 



<j>*r>. 



=A6 2 + 



(19) 



i3 



f=A9 2 + 2^(Vj 2 + v 2 l + v 3 *j +2v 2 2 + 2v 3 l + 2 Vj 2 ,) (19') 

Dans le eas d'un fluide incomnressible, ©» 0, done seul intervient 
le coefficient^ . Le travail de deformation se reduit alors a df . 
La fonction * 4tant positive, on conclut que 1* coefficient de 
viscosity jA est toujours positif. 
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ANNEXE : Expression des composantes du tenseur des contraintes 
et de 1' equation de NAVIER-STOKES en coordonnees cy - 
lindriques et spheriques . 



I) Coordonnees cylindriques 



r, 9 , z 

jites 




Composantes du tenseur des contraintes 

<T = - ■ -» r 

rr 



-p + 2/1 (l/r ^ + 
zz * / o z 



e 



) 



zr sot d z 

L* equation de continuity s'ecrit s 



■r-> 



= 0 



d r r r d 6 * z 

Liquation de NAVIER-STOKES (on neglige les forces de volume) 



^u 



^1 



■+u 



u„ ~iu ~S»U 



© - 



r 

2. 



e ^ r 



^ 2 u 



r , I 



2 iu © 



^6 z ^ 5 



r ^r r 

2V 



»z 



I 



I Vu Vu 



^6 



^u. 



V 



© . I "6 . - 
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/sx z . .. z , 6 z z _ I dp , 

% Vu _ "^ 2 u S 2 u _ ^ 

2 ) Coordonnees spherigues : r , l£ , © 

— t u 

r© / r ^© d r r ' 
W = /* ( rsinfe -^#" + "-~ "Te r 5 



<r - a ( Hi + 1 



Equations de continuity s 



" bUr 1 ^"e I ^ u f 2u r u e cotg 
TT~ + "r" ' rsinB + r + r 
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Equations de NAVIER-STOKES i 

u!T+u 



> u u„ "5u_ u. B -2.. .2 



® ^1 rTr ~r ^8 rsin© r 

I , Si ** 2 <™r> , I *V I H 



u 



r sin' 

, cotg6 ^ U r 2 ^ U e 2 _f^r_ 2cotg " 

<>6 ~7~ T$ r^sinB >f ' 2 " 2 • 

r r _ 



4 6 



r 5 

®TT +,J rlr 1 r 



u m c>u 


u u 
re 


COtgS 


rsin© £f 


r 


r 


V(ru ) j 

© . I © , I 




^r 2 r 2 ^6 2 




2 cos© 5u £ , 2 


K 


u e 7 


r 2 sin 2 G ^£ r 2 " 


"So - 


r 2 Bin^B j 






u © u <£ cot ^^" 


1 rsin<© if 1 


r 


r 


T, ^ 2 (ru f ) z 




c) 2 

i 1 


+ V [ r "br 2 ' r ? 


^e 2 


r 2 sin 2 & It 



, cotge ^ u f , 2 ^ u r , 2cos^ H u if 
r 2 "XI7e W r 2 sin 2 6 " r 2 sin 2 9 
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CHAPITRE V 
V ECOULEKENTS LAMINAIRES ET TURBULENTS 



Introduction : Nombre de REYNOLDS : 

II existe divers regimes d ' 6coulement d'un fluide visqueux. Inexpe- 
rience de REYNOLDS permet de mettre en evidence ces regimes. En i-n- 
troduisant un colorant dans l'axe du tube ou circule un courant 
d'eau, on peut observer les phenomenes suivants : 

-a- pour de faibles vitesses de courant d'eau, le filet de colorant 
traverse tout le tube en se maintenant dans l'axe. 

-b- pour des vitesses plus ^levees, apres avoir accompli un parcours 
plus ou moins long, le colorant se melange brusquement a I'eau. 
-c- pour de grandes vitesses, le filet se melange imiaeu iatement a 
l'eau. 



Dans le premier cas, l'ecoulement est dit laminaire (ou de POISEUILLE): 
les filets fluides glissent les uns sur les autres sans se melanger. 
Dans le troisieme cas, le liquide est anime de mouvements tourbil- 
lonnaires qui produisent un brassage de toute la masse. C'est le re- 
gime turbulent (ou hydraulique) . 

Le deuxieme cas correspond k un regime de transition. REY1I0LDS etu- 
dia le premier ces types d'ecoulement en faisant varier le diametre 
de la conduite, la temperature etc... et montra que le parametre 
permettant de determiner si l'ecoulement est laminaire ou turbulent 
est un nombre sans dimensions appele nombre de REYNOLDS defini par : 



D = diametre de la conduite, v = vitesse,^= masse volumique et 




a 




(I) 
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la viscosite dynamique du fluide. 
On admet en general que si : 

^i<2000 le regime est laminaire 

^>5000 le regime est turbulent 
Ces valeurs limites ne sont d'ailleurs pas tres bien d^finies et il 
faut les considerer comme ordre de grandeur. 

On peut donner une signification physique au nombre de REYNOLDS en 
exprimant le rapport des forces de viscosite" et des forces d'inertie, 
La force d'inertie par unite de masse est : 

Ix dV V__ dV = VdV 



dt V dt dx 



La force de viscosite (voir equations de NAVIER- STOKES (10)) est 

5 M d 2 V 
proportionnelle a 

f dx^ 



Le rapport de ces deux forces est : 

VdV x I 
dx d 2 V 

a la dimension de -j- — 

M d 2 V >» V 

-~- ■ ■■ a la dimension de — — — 

e ~ e i> 

Le rapport de ces deux grandeurs est un nombre sans dimensions for- 
mes par le produit VLp analogue au nombre de REYNOLDS. 

>» 

Le nombre de REYNOLDS est done une mesure de 1' importance des for- 
ces d'inertie par rapport aux forces de viscosity. 

Du point de vue pratique, il est bon d' avoir des ordres de grandeur 
des differents termes qui interviennent dans la definition du nom- 
bre de REYNOLDS : 

Prenons par exemple le cas de l'eau k 20°C ! 

^= 10 3 kg/m 3 , / M= I0~ 3 P1 
Le regime sera laminaire pour(JJ^2000 soit Dv<2 I0~ 3 m 2 /s 
Pour un diametre de conduite D de 0,1m, on doit avoir 
v<2 I?) 2 m/s 



© [G.FANTOZZI], [1975], INSA de Lyon, tous droits reserves. 



- Q 9 - 



Industriellement, les vitesses d ' ecoulements des fluides spn,t de . 



l'ordre du m/s de sorte que les irl^mentSVauront un ecoulement lami- 
naire. Par contre pour des tubes capillaires (D = 10~"Sa), les : 
mes seront en general larainaires (v 6 20 m/s). 



V-I Bcoulements lamlnaires : 

Lorsque les forces d'inertie sont nulles ou negligeables devant les 
forces de viscosite c'est a dire lorsque le nombre de REYNOLDS est 
faible, on peut alors trouver des solutions simples des Equations 
de NAVIEK-STOKES dans le cas d 1 ecoulements permanents de fluides 
incompressible s. Ces ecoulements sont appeles laminaires. 
De faibles nombres de REYNOLDS sont observes dans les cas particu- 
liers suivant : 

- fluide de grande viscosite cinematique 

- vitesses d 1 ecoulement faibles 

- mouvement de fluide de dimensions requites 

V-I I) Ecoulement de P0ISEU1LLE : 

II s'agit de 1 ' ecoulement laminaire d'un liquide visqueux dans un 

tube cylindrique de diametre D. Nous considererons le cas d'un 6- 

coulement permanent d'un fluide incompressible. 

Nous negligerons d'abord la pesanteur. 

^ a _^ Les trajectoires sont paralleles a 

a f \ 

l'axe du tube Oz. 

Les composantes de la Vitesse "^"sont 

done egales k : u = 0 

v = 0 

w 4 o 

Liquation de continuity (^q.II 6) se reduit alors a. : 

divT= 0 soit = o 

Par consequent, w ne depend que de x et y. Par raison de symetrie, 
w ne depend que de la distance v* k l'axe du tube. II est done plus 
commode d'utiliser les coordonnees cylindriques (r,6,*). 
Les Equations de NAVIER-STOKES se reduisent aux equations (voir 
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Annexe chap. Ill) 



(2_ = o) 



0 
0 



r 

- -i- sl. +?Aw = o 

Des deux premieres equations, on en d£duit que p est independant 
de r et de©et ne depend done que de z : p = p(z) 
2 2 w , I 2 2 w + ^ 2 w I ^w 
r" "^© 2 ~dr ? 

w n'est fonction que de r 



Vw 

T + 



So 



at 



I ow _ _I_ d / dw s 
r dr r dr 1 dr ' 

dw 



(r 



) 



k 



Jr ' dr 

La pression p ne dependant que de z et w n'6tant fonction que de r, 
la relation precedente ne peut etre v^rifiee que si chaque membre 
est 6gal a une constante k. 
Integrons cette relation. 

2 



d , dw n kr 

rr- (r -ar-> = + -yr 



dw 



kr 



+ a 



dw 



kr 



a 

r 



kr 



+ a Logr + b 



a et b sont des constantes d 1 integration qui peuvent etre determi- 
ners par les conditions aux limites. 

La vitesse ne peut etre infinie sur l'axe (r=0) done a = 0 et en 

outre, la vitesse s'annule sur les parois du tube (le fluide adhere 

D 



aux parois) : w=0 pour r 



r„ = 



d'ou 



kr! 



Finalement, on obtient pour w 1' expression suivante 

(2) 



k ,2 2, 



k<0 



La vitesse est maximale sur l'axe du tube (r=0) s 



max 



JELL 
4^ 
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w - "max (I " 



-) 



(2') 



La loi de distribution des vitesses suivant un diametre est para- 
bolique (suivant la section, il s'agit d'un parabololde de revolu- 
tion ayant son sommet sur l'axe du tube). 

Remarque : Les calculs precedents 
restent valables si on tient comp- 
te de la pesanteur, a condition de 
remplacer la pression p par la 
pression etoilee p = p + ^gh 
On peut calculer le debit en volume du tube cylindrique. On a 



=fv*dlf =1 °w 21T rdr =' - 




i 



% - - 



2/ M 




( 3 ) 



q = w 

max 



HI* 



La vitesse moyenne est la Vitesse qu'auraient tous les filets flui- 



des avec l€ 


; meme 


d<5bit. 


Elle 


est done ( 




"w = 


-g— = 




w 

max 
2 



(4) 



Perte de charge dans un tube cylindrique : 
Nous avons vu que ■ jjS- = k k<0 

soit p « kz + c , c constante d ' integration 
La pression decrolt done lineairement le long du tube, cette pres- 
sion etant constante dans une section. 



k - dp - 
dz 



L 



(5) 



L etant la longueur du tube et Pj - P 2 = Ap la difference de pres- 
sion entre les sections d 'entree et de sortie. 
Le debit volumique exprime en fonction de 4 p s'ecrit : 
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Tf r 4 a 



(6) 



C'est la loi de HAGEN (1839) et POISEUILLE (1840). 
La chute de pression lineaire s'ecrit done d'apres les relations 
(4) et (6) : Ap _ B^q ^ 



La perte de charge lineaire d'un ^coulement lamfnatpe est par con- 
sequent egale as T P I~ P 2 



(nbre sans dimensions) 



|La charge totale est egale a , u , p , V 2 

Par consequent, si h et V sont constants, la perte de charge le long 
du tube cylindrique est due a la variation.de la pression.~j 



CgD 2 



(7) 



En mettant en Evidence 1 ' energie cinetique de l'unite^ de poids cor- 
respondant a la vitesse de debit w* et le diametre B de la conduit*, 

on peut £crire „ 

J= r^4. (8) ' 

f etant un coefficient sans dimensions appeli coefficient de perte 
de charge lineaire . 

Dans le cas de I'ecoulement de POISEUILLE t 



a 



p wD 



64 



(9) 



etant le nombre de REYNOLDS de I'ecoulement 



Pour caleuler la perte de charge d'un £coulement laminaire le long 
d'un tube cylindrique de diametre D, il faut connaltre w" ou le d6- 
bit et (c'est a dire la temperature et la pression). 

Frottement sur la paroi j 

Le tenseur des contraintes (cf. annexe chap. IV) en coordonnees cy- 
lindriques n'a qu'une seule composante s 



zr 



T U ^ w 
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V-I 2) Ecoulement entre deux plaques planes paralleles : 

Consid^rons deux plaques planes et paralleles (immobiles) , a. la 
distance h = 2e l'une de 1' autre, d'envergure infinie et au re- 
pos. Nous consid^rerons une longueur unite\ L'ecoulement est per- 



manent . 
~7t — 



2e 



L'^quation de continuity (u est la 
seule composante de la vitesse) con- 
duit a la relation : 

=0 u ne depend que de y 
et z. 



Par raison de symetrie, u ne depend pas de z done u n'est fonction 
que de y. 

Liquation de NAVIER-STOKES se r<5duit a : 

r - du * 



du d u c^u Su Su _ 



Ainsi la pression p ne depend que de x et on a liquation suivante 



d 2 u 



ay' 



(13) 



Par analogie avec le cas precedent, = k = Cte 



En integrant 



u = 



+ ay + b 



a et b £tant des constantes d 1 integration 

La vitesse est nulle pour y=eety=-e : par consequent a=0. 

La condition u=0 pour y = e permet de determiner b : 

2 2 
1— + b = 0 soit b = - JE 1— 

La loi de distribution des vitesses est de la forme : 



u = 



k , 2 2* 
(y -e ) 



k<£0 



(14) 



Cette loi est encore parabolique. 
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La Vitesse est laximale au centre : 

2 



ke 



max 



2 / 



h 



Le debit en volume est donne par 1'integrale 



4e 



J-e 



2 

La Vitesse moyenne u = -j 



max 2e 

Les caleuls da t)erte de charge et de frottement se font de ] p rtSme 
facon que precedernment. 

Remarque : L'une des parois peut §tre en mouvement uni orne. Lcb 
conditions aux limites sont alors moditieee et on a alors l'i'coulo- 
ment de GO" ET'fE generalise. 

V-I 3) Notion de couche limite : 

Si on considere une couche de particules tree petites en cot tact 
avec une paroi (qui presente toujours des asnerites) , on no-'it ad- 
mettre une adherence conrolete : la vitesse d ' ecoulevnen t est okc 
nulle pour la coucie dr fluide en contact avec la paroi. 
Quand on s'eloi^ne oe la paroi, la vitesso d * ■' coulenent croft ce 
zero jusqu'a une certaine valeur V qui es 1 , In vitesse d''coulenent 
au loin de la oaroi. Cet accroisserient est d'abord rapide o is e 
plus en plus lent. 

,(\ ( 




fluide litre (la vi' ef.se 
unifor-.e) 



y e i t 



couche liwite 

parti 

La couche de fluide en contact avec la paroi pour laquelle la vites 
se passe de 0 a 0,99 V est amaelee couche Unite dynamique . Au-de- 
1\ de la couche limite, la vitesse est pratiquement constante : 
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c'est la zone de fluide libre ou les effets de visoosite sont tres 



faibles. 



En general, la couche limite dynamique s'accrolt dans le sens de 



l'ecoulement et a partir de l'extr^mite arminV' de la paroi, a cause 
de la transmission progressive des effets de freinage de la paroi. 




Effet de la visoosite : Considerons l'ecoulement d'un fluide de vis- 
oosite variable autour d'un obstacle fixe. 



' f f ' ' /V /x celle qui regne loin de 1' obstacle. 

Si le fluide estvvisqueux et le mouvement tres lent (tres faible 
nombre de REYNOLDS), le profil des vitesses est indique' par la cour- 
be 2 (il y a adherence a. la paroi). 

La r^alite se situe entre les deux extremes (courbe 3 pour les faibles 
valeurs de(Jfe, courbe 4 pour les fortes valeurs) . 

Cas d'un ecoulement permanent dans une conduite cylindrique : 

Considerons un fluide penetrant dans une conduite cylindrique avec 
une Vitesse uniforme ~Tl < 
Quand le fluide s'gcoule dans la canalisation, l'effet de la paroi 
se manifesto : la couche limite s'^paissit jusqu'a atteindre le cen- 
tre de la canalisation a. l'abscisse x . 
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Si le fluide est parfait, le pro- 
fil des vitesses le long d'une nor- 
male en M a la paroi est donn£ par 
la courbe I, la valeur Ves £tant 



- 10? - 



A rartir de I'abscicse Xg, la zone de flui e libre est reduite a 
nennt i le profil des vitesses a alors une allure parabolique si 
1 * ecoulement est a^sez lent (ecoulement laminaire) ou plus aplatie 
si 1 ' ecoulerent est turbulent. 

A f -./tir d<=> x a , le profil des vitesses n'evolue pas : le regime est 
etabli . Cette lonrueur x^ defend du nombre de REYNOLDS. Pour les e- 
conle ients mdustriels (turbulents) , x^ est de I'ordre de grandeur 
de 'ilusieurs dizaines de dia'netrec. 
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V.2.- ECOIIEMENTS TURBULENTS 

La plupart des ecoulements rencontres industriellement etant turbulents, 
il est done important du point de vue technique d'etudier ce type d' ecoulement. 

V.2.I.- CARACTERISTIQUES DE E' ECOULEMENT TURBULENT 

Dans un ecoulement turbulent , la vitesse d"ecoulement en un point M 
subit des fluctuations continuelles de petite anplitude en direction et en 
intensite autour d'une valeur moyenne independante du tenps(si les conditions 
f aux lindtes sont constantes 3 on a 

A done un ecoulement permanent en 

moyenne) 

Remarque : dans le cas de 1' ecoulement 
laraLnaire, 1' ecoulement est rigoureu- 
sement permanent (les fluctuations 
sont indecelables S la:precision des 
mesures pres) . 

a) DefinitiaiS : On appelle fluctua- 
tion u' de u(t), la variable egale 
a la difference entre chaque valeur 
individuelle et la moyenne 




u' = u(t) - u (15) 



La variance e st la moyenne du carre de la fluctuation u 1 et 1'ecart type la 
racine carree 0"(u) de la variance : 



CT 2 (u) = u' 2 = [u(t) - u]' 
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La covariance de deux variables u et v est la moyenne du produit de leurs 
fluctuations : 

u' v ( 

, La ncyenne d'une fluctuation est nulle. En effet, la moyenne de la vitesse 
prise entre les temps t. et t ? - t. + T est egale a : 



fc 2 



u(t) dt = u (T, t 1 ) (16) 



II existe une duree seuil de la medi a tion T telle que la moyenne u soit 
constante, quel" que soit t^ 

T ? fc s 

Mnai, on peut ecrire : 

■ - 1 I u dt = i { (u + u') dt = ^ / u dt + | u' dt 



u 



t 1 > t t t 1 /t 1 



soit u = u + u' 

d'ou u>~ = 0 (17) 

« La moyenne d'un produit differe en general du produit des moyennes : 
u v ^ u v 



en effet u v = (u + u' ) (v + v 1 ) = u v + u' v + u v' + u' v 1 
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La moyenne d'une somme est egale a la somme des moyenne et la moyenne du pro-i 
duit d'une constante par une variable est egale produit de cette constan- 
te par la moyenne de la variable. 



On a done : u v = u v u'v = v u 1 =0 

u v' = u v' = 0 



soit u v = u v + u' v' (18) 
Si on egale v a u, on obtient 

T = u 2 + ^ 

soit u' 2 = cr 2 (u) = u 2 " ~ u 2 (19) 
• La moyenne de la dgrivee partielle est egale a la d£ii*#e de la moyenne 

HL*! (20) 

3x 3x v ' 



Eneffet : {- (i J * u «) - J J * g- dt 



fc l et ^2 ® tan * independant de x. 



« La moyenne de la derivee d'une variable par rapport au temps sera supposSe 
nulle. 
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La difference u(t 2 ) - u (t^) est faible et si T est assez grand. 



(21) 



b - Contraintes dans le eas d'ecoulements turbulents 

Nous ne considererons que le cas des fluides incompressibles dont le 
mouvement moyen est permanent. 

Nous allons determiner le tenseur des contraintes dans le cas d'un ecoule- 
ment turbulent 5 partir du tenseur etabli au chapitre IV (relation IV. 8) . 
Aux contraintes moyennes en un point M 5 il faudra aj outer des terroes supple- 
msntaires dus aux fluctuations. 

Nous allons montrer que le tenseur des contraintes a la forme suivante : 



«£. =-P + 2 y || - p u' 2 



'11 



a 22 = - p + 2 " w " pV 



T - , „ aw -.2 

J 33 = " p 2 v a? ~ pw 



T ,au , 3V N —7—7 

"12 - » <» * i5> " » 



V 



T ,av r awv — s ? 

°23 = " ( H + 3? } p V W 



T _ ,3U j 3W, — 5— r 



13 V 3Z 3*' 
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ou sous forme tensorielle 



p 5 ij + 2v v ij " p u 'i u 'j 



(22) 



3u. 



3u 

2V « = 3x7 + 3xJ 



vitesse de deformation 



Du point de vue dynamique, un ecoulement turbulent se comporte comrae un 

(Scoulement permanent ordinaire en imaginant un mouvement fictif qui a les 

earaeteristiques moyennes de 1' ecoulement u, v, w, a condition d'aj outer 

aux contraintes moyennes des termes supplenentaires dus aux fluctuations 

de la forme p u*. u*. 

j 

Etablissons maintenant les relations (22). Les composantes de 1' acceleration 
dans le cas d'un fluide incompressible s'ecrivent de la facon suivante : 



du. 
l 

dt 



3u^ 3u^ 
3t + U j 3x7 



1' equation de continuity s'ecrit : div V = 0 soit 



3x^ 3X 3y 3z 3Xj 



3u. 3U. 3U. 3U. 3(u.U. ) 

Oi peut done ecrire y ± = ^ + Uj ^ + Uj_ ^ = ^ + 



soit par exemple 



3u 3u 3uv + 3uw 
Y x ~ 3t 3x 3y 3x 
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Le mouvement moyen etant permanent, on peut ecrire 



_3u + 3v + 3w _ q 
3x ay 8g 



et 



3U i U j 
8X j 



3U i U j 
3X. 



d'apres la relation (20) 



u.u. = u. u. + u'.u'. "d'apres (18) 



au. u. 
3xi 



soit 



3u. u. 3u'. u' . 
1 ,T + 11 J 



ax. 



3X- 



ainsi 



— _i>u au" . au w , au*^, au*v' L au'w' 

Y - — + — — + — — + + — — + 

"x 3ot 3y 3Z 3X 3y az 



au. u. 

Le ler terms — - — i 



represente la composante r x de 1' acceleration du 
mouvement moyen J (d'apres les relations precedentes) . 
Examinons maintenant les equations dynamiques d'un fluide visqueux 



1x7 + p x i - p It - 



Prenons les moyennes dans le temps : on obtient : 

3u.u. 



3o. . 

axj 



3u'.u'. 

pX = pv. = p -ri-i + 0 — i 

4. l 3X4 p a x - 



J 



r i 



,-i — 3° ii 



au'.u'. 



3x,- 
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Cette equation peut s'ecrire sous la forme suivante (identique a celle 
d'un fluide visqueux) . 



avec 



3aT. 



<j. .= a. . - o u' .u. 

1J 1J 9 i J 



(24) Equations dynamiques du mouve- 
nEntjnioyen d'un fluide visqueux 

turbulent 

(22) 



T 5 ij 



+ 2p v. 



3U. 3U. 
D{., + 2u (-r-i + — *4 

F ij p v 3x. sx^ 



On retrouve bien 1' expression (22) que nous avons^au debut du paragraphe . 
Les equations dynamiques du mouvemsnt raoyen d'un fluide visqueux turbulent 
inconpressible s'ecrivent (voir paragraphe IV.2.2.) 



i = x i"F Sr + ^i-<— (2 5) 



j 3u' 2 + 3u'v' + au'w ' J 
I 3x 3y 3Z | 



Ces equations sont appelees equations de REYNOLDS 

Aux equations dynamiques, il faut aj outer 1' equation de continuite soit 
pour un fluide inconpressible : 

divV = 0 = |H + H + |« = o 
3x 3y 3z 
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Soit en considerant les valeurs moyennes 




car 



u 



i = 



w 



' = 0 



a chaque instant on a d'ailleurs : 



3u' 
3x 




V.2.2.- DISTRIBUTION PES VITESSES 

Les composantes suppleroe ntair es du tenseur des contraintes dues i la 
turbulence sont de la forme pu' 2 pour les tensions normales et pu'v' pour 
les tensions tangentielles , Ces tensions sont ggneralemsnt beaucoup plus 
grandes que les forces de viscosity si bien qu'on peut negligsr ces dernieres 
dans les equations dynairiques. 

Pour resoudre ces equations, il faut expriner u' , v' etc... a" partir de 
certaines hypotheses. Considerons le raisonnement de PRANDTL. Soit un ecoule- 
ment plan ou la vitesse moyenne, parallele a l'axe des x, ne varie qu'avec 

y- 

PRANDT suppose que la fluctuation u' est 
proportionnelle a la pente du profil des 
vitesses : 




(26) 



v' etant donne par une expression analogue 
la longueur 1 est appelee longueur de 
melange et caracterise la turbulence. 
Ainsi la tension tangentielle de 
turbulence est egale a : 




did 




M- 
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Si on compare cette formule avec la formule de NEWTON x = v ^ 
qui donne la conposante visqueuse, on peut considerer que le coefficient 

Pi 2 | | | (27 

appele coefficient de viscosite dynamique de turbulence joue le meme role 
que p. mais depend de y contrairement a (en general oe coefficient est 
» vi ). 

Par 1 des considerations de similitide, VON KAJMAN a ete amsne a" poser : 



k Stant une constante universelle. 

Dans certains cas sirples, on peut effectuer des hypotheses qui perrcettent 
de rSsoudre les Equations. Alhsi au voisinage d'une paroi solide (plaque plane, 
tuyau, etc.), on peut considerer que 1 est ppoportionnel a" la distance y 
a* la paroi : 




(28) 



1 = Cy 



soit t t = p C 2 y 2 (|i) 



Si on suppose que x^ est constant et egal a" x Q sa valeur sur la paroi ( 
(y=o) on a, 




soit 




(29) 
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Cette formule represente assez bien la repartition des vitesses moyennes 
observees experlirentalement (sauf au voisinage de la parol : couche limite 
turbulente) . 



Intensite de turbulence 



cCty Composante. u de la vitesse suivant l'axe 0 X varie avec le tenps au 
tour de sa valeur moyenne u. 



4^ 



t 



L' importance de la turbulence est 
caracterisee par 1' expression 



MS 

V 



(30) 



qu'on appelle i ntensite de turbulence 
suivant Ox (V etant la vitesse 
rooyenne de l'ecoulement) . 

L' intensite globale de la turbulence 
sera definie par la relation : 



^(u' + v' + w' ) 



Bnergje du mouveirent turbulent : 

A un instant donne, 1'energie irecanique totale d'un volume unite de fluide 
est egale 5 : 

E = p +j> gi + j 1 y- 

Mais la pression et la vitesse subissent des fluctuations 

P = P + P' 
V = V + V 
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Calculous la valeur moyenne de E : d'aprea la relation (19) 

5 = Pr +p + p^i + p-~ = E 1 + E"' (3D 

aveo E* = p — = (u' 2 + V 2 + w' 2 ) 

representant l'energie totale du irouvement rooyen 
W representant l'energie cinetique de turbulance 

Lorsque la turbulence est etablie dans un ecoulement, u 1 , v^, w' ne 
varient plus et par consequent E"' : 1' ecoulement est alors 3tationnalre . 

Remarque : dans un ecoulement reel (laminaire ou turbulent), il y a toujours 
diasipation d'Snergie mecanique a cause du frottement interne. Dana le cas 
du mouvement laminaire, cette dissipation se fait Si l'Schelle mol&culaire 
entre les couches de fluide. Dans le cas de l'ScouLement turbulent, en plus 
de ce processus moleculaire, s^aAOwte on processus S une Schelle plus 
importante : il y a echangs entre particules fluides de volume plus ou moins 
important. 

Ces particules qui passent d'une coucbe fluide a une autre mettent en jeu 
des forces bien plus importantes ; ainsi les phenomenes d'origine purement 
visqueuse sont masques par les phenomSnes d'origine turbulente. 



© [G.FANTOZZI], [1975], INSA de Lyon, tous droits reserves. 



- 119 - 



V.2.3-- PROFIL PES VITESSES ET PERTES EE CHARGE 
1 - Charge dans une section 

Nous avons defini la charge totale par 1 'expression : 

^ Bh+ {S + -n (32) 

La charge tbtale represente, en hauteur de fluide, l'energie mecanique totale 
de l'unite de poids du fluide (P = mg = 1) . 

L' expression p— + p + pgh represente, en unites ''de pression, l'energie 
mecanique totale de l'unite de volume du fluide. 

Quand on considere le mouvement permanent d'un fluide parfait incompressible, 
le long d'une trajectoire, (il satisfait a la formule de EERNOULLI) . 

V 2 v 2 
P — + P 1 + P&i ± = P + P 2 + Pgh 2 

ou h_ = hm : la charge totale est constante le long d'une 

x l 2 

trajectoire : ceci traduit la conservation de l'energie mecanique de la par- 
ticule entre 2 positions 1 et 2. Si le fluide est visqueux, on ecrira pour 
une particule suivant sa trajectoire : (equation IV. 1 2.). 

dr etant le travail mecanique perdu par le fluide (transforme en chaleur 

par frottement visqueux) par unite de volume du fluide. 

Integrons cette relation entre les positions 1 et 2. 

2 2 

V * V ? * 

P — + P, = P ~ + P 3 + t (33) 

I 1 1 2 

t = travail effectue par les forces de frottement par unite de volume 
(perte d'energie mecanique) . 
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Dans le cas de fluides liquides, on exprime la relation de BEHNOUILLI en 
hauteur de liquide : 



ll + h*- Y JL*JL 

2g Pg 2g pg 



Pg 



la difference Ahp = h^, - hp = 



T 

Pg 



represente la perte de charge entre les deux positions (1) et (2) de la 
trajectoire de la particule. 

On peut definir deux types de pertes de charge : - les pertes de charges 
regulieres ou reparties : il s'agit des pertes de charts dans des conduites 
longues (l'gcoulement gtant soit laminaire (relations V 7, 8 et §) soit 
turbulent). 

-les pertes de charges singulieres : elles sont dues I des singularites ou 

accidents existant dans les conduites (elles se calculent 3. l'aide du theo- 
remev: voir chapitres II et III). 

Nous avons de"fini pricedeitment la perte de charge linSaire par la quantite 
J : (relation V.7). 



h Tl - hm 
L 



(35) : nontore sans dimensions 



et un coefficient de perte de charge lineaire v f' defini par (voir relation 
V.8). 



vi^ u 1 
I 2g D 



(36) 



(Ce coefficient fait apparaitre 1'energie cinetique de 1' unite de poids et 
:Le diamltre D de la conduite) 
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2 - CALCUL IE LA. PEKIE DE CHARGE 

a - position du problems 

Si 1'ecoulemsnt s'effectue dans des conduites circulaires lisses, la 
perte de charge entre deux sections distances de L peut s'exprluer par la 
relation- suivante : 

A *= " pj = f (V, D, L, Pr u) (37) 

En effet, la conduite etant cylindrique, la section est constante et la 
vitesse iroyenne est constante . 

Des considerations d 'analyse dimsnsionnelle ucntrent que l'on a : 
_2 

f p "2" D f ( ~ } 

ap*= p 5j I f my (38} 



SOlt 



* ¥ 

P t - P 2 v L 

- — — 2i — f ^ 



Le coefficient de perte de charge lineaire est egale a. f ((ft) 

^ = f «R> 

Ce coefficient suffit pour definir la perte de charge : il ne depend que de 
<R. 

La determination theorique de f(fl) est facile dans le cas de 1'eeoulenent 
lamLnaire mais pour un ecoulement turbulent, il faut faire appel 

a 1' experience. 
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Dans 1' Industrie, les conduites ne sont pas lisses mais rugueuses. II faut, 
par consequent, tenir eompte de cette rugosite. Aussi, la relation (37) 
sera remplacee par la relation suivante: 

* k 

f 1 - P 2 = f (V, D, L, p, v, rugosite) 

La rugosite est une notion dirflcile a definir : elle peut etre definie par 
plusieurs facteurs. On peut essayer de la caracteriser par un parametre moyen 
coimie la hauteur moyenne £. dee aspirites : 

Ainsi p 1 - p 2 = f (V, D, L, j> ,V ,£) 
tia raisonnertfint analogue au precedent permet d'ecrire : 

bf'fir^**^ (39) 

Le coefficient de perte de charge y = f (IB., ^) est^du nonbre de REYNOLDS et 
de la rugosite 

Nous pourrons determiner la relation qui existe entre le coefficient 

et le coefficient de frottement interieur C f defini par la relation V.ll : 

'p ■ Ml" 2 

Tp etant la force de frottement par unite d'aire sur la paroi 
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Considerons le fluids contenu entre 2 sections 1 ep 2 et appliquons le theo- 
rems d r EULER a la surface de reference constituee par les sections 1 et 2 
et les parois de la conduite. 



I F = J V6^ d<| = debit de masse 

systeme forces appHquees au fluide contenu dans S 



Le fflouvement etant permanent 



= 0 




Par consequent on a : 
IDL 



<pf- 



4 > s 



Remarque i pour un fluide pes ant, 
on doit considerer les p* au lieu 
des pressions p pour un fluide non 
pesant. 
Aijnsi : 



"f 2 



nDL 
S 



s = 



IID 



Soit : 



= tc f p 



5 d ' oil 



* = 4 c. 



(40) 



-b- experiences de NIKURADSE 



NIKURADSE a etudie les variations de f en fonction de J? et f , pour des 
conduites avec rugosite homogene (obtenue a l'aide de grains de sable calibres) 
II a obtenu les resultats indiques par la figure suivante : 
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Bout une valeur donnee de la rugosite e/p, les variations de $ aveo fff peuvent 
etre representees par 3 portions de droite (J) , (?) gk Q) qui se racoordent 
entre elles par des zones intermediaires . Quand ^ varie, la position de la 
droite @ varie tout en restant horizontale. 
Pour une valeur donnee de p on peut definir cinq regions : 

- les rggions I ; II et III ou ifi ne depend que de <R : on dit que la conduite 
est hydrauliqueiiBnt lisse. 

- la rSgLon V : i|i ne depend plus que ■—- : le frottement est rugueux. 

- la region IV est une region de transition : le frottenent y est seni-rugueux. 

Dans la rigion I, l'ecoulement est laninaire et obeit a la loi ie HACSH- 
PoiseUH.t£ (1) (valable jusqu'a Ri 2000). 



64 



(41) 
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La region II est une region de transition mal definie qui correspond a 

1' apparition de l'ecoulement turbulent. 

La region III correspond a la droite de BLASIUS (2) 

+ = = (100«?)" 1A (42) 
[ 

valable jusqu'a. (R = 10^ 

Si (R est > 1C? , les points experimentaux se situent legerement au-dessus 
de ia droite de BLASI US (courbe (fj)> : on a alors la courbe de RAHMAN : 



Dans la rSgion V, * est donnee. par la formule de KAREMfiN-PRANDTL 

_1/2 = 2 log 10 (3,71 1) (44) 

Ainsi ie comportenent d'une conduite est assez eomplexe : il depend a la 
fois de R (en particulier du debit) et ^ . 

Cas des conduites industrielles 

Dans les conduites industrielles, les aspirites ne sont pas regulieres : 
la rugosite est heterogene. Toutefois, les experiences de C0LEBR00K ont 
montre qu'on pouvait generaliser les resultats obtenus dans le cas de conduites 
dont la rugosite est homogene , En particulier pour d? < CR^ la conduite est 
hydrauliquement lisse alors que fft > ff^ la conduite est hydrauliquement 
rugueuse. La relation (44) permet alors de deflnir une rugosite uniforms 
equivalente £. . 

Le tableau suivant donneles valeurs de c pour un certain noirbre de parois 
usuelles . 



1/2 



2 log. 



(R * 



1/2 



10 2,51 
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Matgri.au 


£ (mm) 


tube etire cuivre 


0,001 


acier neuf 


0,045 


acier rouillS 


0,2 


acier biturae 


0,015 


fer galvanise" 


0,15 


fonte 


0,25 


fonte goudronnee 


0,1 


b§ton 


0,3 a 3 


planohe en bois 


0,2 1 0,9 



Dana la pratique, il faut tenir conpte du fait que la rugositfi d'une 
conduite varie avec le temps. ^— - 
Si on porte sur un graphique - 2 log r- en ordonnee et 2 log -fi- 
eri abscisse, les resultats obtenus avec les conduitro industrielles se 
situent sur une courbe unique representee par la formule enpirique de 
COIEBROOK 

4 - 2 log (£) = 1,14 - 2 log (l ♦ 9,32 J^) (45) 



ou 



Cette relation (45) est la formule universelle des pertes de charge dans les 
conduites indsutrielles 
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f Mr 




On peut constater que las conduites industrielles ont un conportement 
analogue a celles d« NIKURADSE pour las petites et les grandes valeurs de 

. Dans la pratique, le calcul de V se fait 8 I'aide d'abaques (par 
exenples le diagramre de MDODY qui donne en fonction de (R pour diffe- 
rentes valeurs de Ce diagramrre est identique a celui de NIKURADSE 

sauf dans les zones de transition. 
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3 - Profil des vltesses 



Dans le cas tie l'ecoulement laminaire, la loi de repartition des 
vitesses est parabolique. Exanrinons maintenant le oas de l'ecoulement tur- 
bulent. 

La vitesse u s'exprirae par une relation de la forme 

u = f (f pJ D, p,p ,e , y) (46) 

t Stant la force de frottement par unite de surface sur la paroi. 




u^, a les dimensipns d'une vitesse et est appelee vitesse de frottement 
Nous pouvons montrer que (conditions de similitude) 




(48) 
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Conduites lisses : la eonduite est hycirauliquemeftt lisse si : 



V 



< 5 



Le profil des vitesses peut etse divise en deux rSgions : 

• d- region voisine de la paroi ^0,15^ : la repartition des vitesses 

ne depend pas de D et £ ; on a :° 

Cette region est constitute de plusieurs zones. Centre la paroi; il y a 

un film landnaire, de tres petite epaisseur . On. peut admsttre que la formule 

de NEWTON s 'applique : 

x = v ^ = constante = 
SOit T = y h Za = v % 

p y p y 



u 


yu f 


u f 





(50) 



Dans le film laminaire, la repartitition des vitesses est lineaire tant que 



y u- 



<5 



y u f 



Si — — est^-25» la repartition des vitesses est logarithmique 



5,65 log 



y u„ 



+ 4,9 



(51) 



Bitre les deux zones, il existe une zone de transition. 
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2 - Dans la region eentrale, la repartition des vitesses ne depend pas de la- 
nature de la paroi. L' experience montre que : 



u . . . - u 
max 



= f \ L ) 



Pour - — J 0,2, on verifle assez bien la loi parabolique 
o 




(52) 
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Conduites rugueuses 
eU-p 

Lorsque est ^ 70, la conduite devient rugueuse et le profil des 

vitesses ne depend alors que de e (et plus de R). 

Inexperience conduit a la loi de distribution des vitesses pour la region 
voisine de la paroi : 



= 5,65 log 2. + 8,1 



(53) 



Si 5< -— < 70, le frotterrent est seni-rugueux et la distribution des 
vitesses depend a la fois de la viscosite (R) et de la rugosite e : 

Dans la region centrale, les memes lois que pour les conduites lisses sont 
applicables (relation 52) 
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